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frequency  range  bettec^han  0.01  to  100  Hz  and  responds  to  minimum  temperature 
fluctuations  of  0.015'u.  Of  special  importance  is  the  low  construction  cost 
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I.  INTRODUCTION 


In  stud) ing  the  effects  of  the  turbulent  atmosphere  on 
optical  propagation,  wc  must  measure  the  strength  and  character 
of  the  turbulence  in  order  to  reference  propagation  measurements 
to  the  current  conditions.  Without  such  a measurement  of  con- 
ditions, the  propagation  results  would  only  be  phenomenological, 
and  it  would  be  difficult  to  extend  the  measurements  to  other 
sites.  It  is  not  sufficient  to  measure  gross  meteorological 
parameters  such  as  wind  speed,  wind  direction,  humidity,  baro- 
metric pressure  and  temperature;  although  these  are  important. 
Rather  microtemperature  fluctuations  are  sufficient  to  describe 
the  turbulence  which  affects  optical  propagation. 


Refractive-index  depends  on  atmospheric  pressure,  water 
vapor  and  temperature.  Pressure  fluctuations  are  of  suffi- 
ciently low  amplitude  to  be  totally  ignored.  Water  vapor,  or 
humidity,  fluctuations  are  significant  in  the  radio  wavelengths, 
but  in  the  optical  range  their  effect  over  land  is  typically 
less  than  10%.  So  we  are  left  with  temperature  f luctutations , 
and  that  is  what  we  measure  and  relate  directly  to  refractive- 
index  turbulence.  Refractive-index  is  not  measured  directly, 
except  to  the  extent  that  we  also  propagate  optical  waves 
through  the  atmosphere  and  record  the  effect  of  the  rapidly 
clianging  refractive-index. 


The  basic  descriptor  of  temperature  turbulence  is  C^, 
the  temperature  fluctuation  parameter,  which  is  defined  as  the 


.1, 


constant  of  proportionality  in  tlie  inertial -subrange  form  of 
the  temperature  structure  function  D^Cr) , given  by  the  Kolmo- 
gorov form 


IV(^)  = <CT  -T  )2> 


= a <<r<<L  , 

1 0 0 


(1) 


where  and  = temperatures  at  positions  1 and  2 with  separa- 


tion r,  - inner  scale  and  = outer  scale.  Notice  that 


is  not  defined  from  measurements  of  where  r is  outside  the 


range  indicated;  and  moreover  if  the  exponent  of  r is  not  2/3, 


then  the  turbulence  is  not  Kolmogorov  and  is  not  defined  at 


all.  We  shall  consider  an  instrument  which  will  give  us  mea- 


sures of  as  well  as  functions  such  as  D^(r)  and  the  tempera- 

T 


j --  --  --  .-^1 

ture  fluctuation  temporal  power  spectrum  F_(f). 


Our  basic  requirements  on  the  probe  are  that  it  be  suf- 
ficiently short  to  resolve  small  scale  fluctuations,  and  that  it 
be  of  sufficiently  small  diameter  to  be  able  to  respond  rapidly 
enough  to  the  highest  frequency  fluctuations.  Also  the  sensor 
circuitry  must  respond  to  very  low  frequencies  to  account  for 
fluctuations  of  size  larger  than  the  outer  scale.  Typically  the 
inner  scale  is  a few  millimeters  and  the  outer  scale  is  on  the 
order  of  height  above  ground.  We  can  use  these  scales  to  desig- 
nate a desired  frequency  range  of  v/(2ttL  ) to  v/(2Tri,  ) where 

0 0 

V = wind  speed.  If  we  suppose  v ranges  from  1 to  5 m/sec, 


a ^ 6 mm,  and  = 10  m,  then  the  nominal  frequency  range  is 


,fi^ 

■i 


- t 


10-^  to  10^  Hz.  Our  design  will  have  a somewhat  better  fre- 
quency range. 

The  particular  sensor,  with  its  associated  circuitry,  is 
described  in  the  next  section.  An  overall  view  of  the  unit  is 
provided  as  fig.  1.  In  general  the  probe,  which  is  a short, 
extremely  thin  platinum  wire,  is  mounted  as  one  leg  of  a Wheat 
stone  bridge  circuit.  The  resistance  of  the  probe  is  very 
nearly  proportional  to  temperature;  and  thus  temperature  changes 
are  sensed  as  an  imbalance  voltage  of  the  bridge.  This  imbal- 
ance voltage  is  amplified  so  that  the  overall  sensor  has  a re- 

sponse  in  the  order  of  0.5  v/“K. 

We  shall  follow  the  sensor  description  of  the  next  sec- 
tion  with  instructions  on  how  to  manufacture  probes,  an  outline 
of  the  steps  employed  in  field  use  and  detailed  instructions  on 
how  to  manufacture  a complete  microthermal  sensor  unit.  Then  we 
will  conclude  with  a description  of  the  data  processing  proce- 
dures  used  in  the  current  operation. 
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II.  SYSTEM  DESCRIPTION 


A.  The  Probe  Itself. 


The  microthermal  probe  is  a piece  of  IVollaston-process 
silver-on-platinum  wire  which  has  been  etched  so  as  to  remove 
the  silver.  The  remaining  platinum  acts  as  the  active  element. 
The  Wollaston  process  allows  the  drawing  of  extremely  thin  wires 
by  inserting  a platinum  core  inside  a much  larger  tube  of  silver. 
The  silver  is  of  a diameter  100  times  larger  than  the  core.  The 
composite  is  drawn  in  a conventional  manner  and  both  materials 
are  reduced  in  diameter.  The  resultant  wire  is  roughly  250  ym 
in  diameter,  with  a platinum  core  of  about  2.5  ym.  A short 
section  (about  1 to  1.5  cm  in  length)  of  the  wire  is  soldered 
at  its  ends  to  a pair  of  posts.  Then  about  a 2 mm  length  of  the 
wire  is  immersed  in  nitric  acid  to  remove  the  silver.  The  re- 
sultant resis.tance  of  the  exposed  platinum  is  about  20  n/mm. 

For  a photograph  of  the  probe,  refer  to  Fig.  2. 


Platinum  has  a very  nearly  linear  temperature -resistance 

characteristic  curve,  with  a constant  of  proportionality  of 
about  0.00393  n/‘’K.  More  precisely  we  can  describe  this  curve 

with  the  Calender-van  Dusen  equation*: 

Rs(T) 


Rjo) 


= 1 + a[T+6(l-T/100) (T/lOO)  + B (1 -T/lOO) (T/lOO) M (2) 


where  R^(T)  = sensor  resistance  at  temperature  T in  »C,  and  for 
pure,  annealed,  strain-free  platinum  wire  we  have: 


..'.1 . . 


Figure  2.  Photographic  view  of  probe  end  of  Microthermal 

Sensor.  Cap  has  been  removed  and  protective  sleeve 
pulled  back  to  put  probe  in  its  normal  operating 
configuration . 


a = 0.00393  f^/“K 
6 = 1.49 


and  B = 0.110  for  T<0  and  B = 0.00  for  T>0. 


For  convenience,  values  of  Hq.  (2)  for  the  temperature  range  of 
interest  are  given  in  Table  1. 


Ochs^  uses  sensors  which  are  0.6  Pm  in  diameter,  and 
they  have  a response  time  of  ahc’;t  200  ysec  in  still  air  or  about 
70-80  ysec  in  winds  of  2 to  5 m/sec.  Measurements  of  the  re- 
sponse time  of  a 2.5  ym  wire  were  made  by  Ochs®  as  900  ysec  in 

still  air  and  300  ysec  in  winds  of  about  5 m/sec.  Wyngaard** 
corroborates  these  values  with  measurements  of  the  response  time 
of  a sensor  used  in  the  RADC  experiments,  finding  about  800  ysec 
in  still  air.  The  values  are  for  "dirty"  sensors,  ones  which 
have  been  used  for  a day  or  so.  "Clean"  sensors,  freshly  made, 
may  have  response  times  at  least  one-half  the  aforementioned 
values.  By  response  time  we  mean  the  time  it  takes  the  output 
voltage  of  the  bridge  to  go  63.2%  of  its  final  value  after  the 
sensor  sees  a step  change  in  temperature.  Since  the  response  is 
exponential,  the  sensor  acts  as  its  own  low  pass  filter  with  a 
power  response  of  [1+  (2TrfT)M’'  where  T is  the  response  time. 


B, 


The  Circuit. 


The  sensor  is  mounted  as  one  leg  of  a Wheatstone  bridge 
shown  in  Fig.  3 and  is  indicated  by  RgCT).  The  balance  poten- 


tiometer Rp  is  adjusted  so  that  the  output  voltage  is 


Table  1. 


I^sCT)/!'.  (0)  versus  temperature  T(°C)  for  typical  pure, 
annealea,  strain- free  platinum  resistance  temperature 
From  Rosemount.* 


sensor. 


T(°C) 

ns(T)/R^(0) 

0 

] .0000000 

1 

1 .0039830 

2 

1.0079648 

3 

1.0119455 

4 

1 .0159251 

5 

1.0199035 

6 

1.0238807 

7 

1.0278568 

8 

1.0318316 

9 

1.0358053 

10 

1.0397778 

11 

1.0437492 

12 

1.0477193 

13 

1.0516882 

14 

1.0556561 

15 

1.0596228 

16 

1.0635883 

17 

1.0675527 

18 

1.0715157 

19 

1.0754777 

20 

1.0794385 

21 

1.0833982 

22 

1.0873567 

23 

1.0913139 

24 

1.0952699 

25 

1.0992249 

26 

1.1031787 

27 

1.1071313 

28 

1.1110828 

29 

1.1150330 

30 

1.1189821 

To  convert  a measured  value  of  Rg  at  T to  Ks  at  0°C,  divide 
appropriate  number  above. 


nearly  zero.  (This  is  done  when  tiie  switch  is  in  tlie  DC  posi- 
tion.) The  imbalance  voltage  clue  to  temperature  fluctuations  is 
amplified  by  a series  of  two  operational  amplifiers  shown  in 
Fig.  3.  The  circuit  is  sufficiently  straightforward,  so  we  will 
not  go  into  its  details.  Later,  under  construction  details,  we 
will  provide  photographs  as  well  as  printed  circuit  board  lay- 
outs of  the  circuit. 

For  the  l)ias  voltage  and  resistance  values  given  in 
Table  2,  and  applying  to  Fig.  3,  we  find  that  the  current  drawn 
by  the  probe  is  about  450  pa.  According  to  tests  made  by  Ochs^ , 
this  implies  we  have  a cold-wire  probe  as  long  as  the  wind  ve 
locity  exceeds  about  0.5  m/sec.  In  fact  we  recommend  th.^t  the 
reader  peruse  that  report  for  significant  supplemental  informa- 
tion. 

The  sensor  amplifier  has  a pair  of  capacitors  each 

labeled  C between  first  and  second  stages  to  allow  AC  coupling 
2 

of  the  stages.  The  values  used  provide  a 100  sec  time  constant 
and  thus  the  -3  dB  point  of  the  high  pass  single  pole  filter  is 
at  1.6x10'^  Hz.  This  is  sufficiently  low  to  encompass  scale 
sizes  larger  than  most  values  of  L^  and  for  typical  winds.  As 
long  as  the  circuit  is  enclosed  in  a nearly-airtight  container 
and  the  components  have  low  drift,  there  will  not  be  a signifi- 
cant low  frequency  drift  problem. 

A relatively  stra"' ghtfonvard  analysis  of  the  circuit  in 
Fig.  3 gives  us  this  gain  equation: 

0 


"ode 


A1 
A 2 
Cl 
C2 


C3 

C4 


C? 


>3 

Ra 

Rb 

R1 

R2 

D 


„) 


R .1 


Rp 


Table  2 


Quantity 


R3CT) 


RiM  SW 
CAL  SW 


Cable 

Cable 


1 

1 

1 

7 


1 

1 

1 

2 

1 


3 

1 

7 


1 

1 

7 


1 spool 


1 

1 


Component  list  for  Bridge/.Amplif ier 


Value 


500 

200 


yfd 

yfd 


390  yyfd 
220  yyfd 
10  yyfd 
30.1  Kn 
2.74  Tn 
500  Kfi 
1 K9. 
10  K9 
100  9 

1 m. 

200  9 


as 

as 

as 


req ' d 
'd 


req 


req 


Item 


±5?6 

±5?o 

±n 

±n 

±1% 


±v 


0 

±u 


±Vo 


±1*' 


IRC 

IRC 

IRC 

IRC 

IRC 

IRC 


Resistor 

Resistor 

Resistor 

Resistor 

Resistor 

Resistor 


Analog  Devices  A0208I!  Op.  Amp. 

Analog  Devices  M)504U1  7235H  Op.  top. 
15  WVDC  Mallory  MTV500  DW15  Capacitor 
25WVDC  Sprague  TR-1213,  Capacitor 
30  D207G025DH2 

+5%  Elmenco  Mica  7FA391G031  Capacitor 

" Elmcnco  Mica  7FC221G03  Capacitor 

Elmenco  Mica  6CD100K03  Capacitor 

IRC  Resistor  T055  (MIL  RN55D) 

TO  5 5 

T060  (MIL  R1J6  0D) 

T055 
T055 
TO  5 5 

zi-a  iixv.  T060  QQ.QTon-i 

Amphenol  Potentiometer  994SL201 
Amphenol  Connector  7 Pm  #197  (Panel 

Plug)  , O TPC 

IC  Socket  Cinch  - Jones  8-1  Lb 
Wire  0.1  mil  diameter  platinum, 
silver  coated,  5747  f2/Ft  at  20  C, 
H08,  Sigmund  Cohn  Corp.,  -'t. 

Vernon,  NY  „ . 1 

Subminiature  Switch,  Switchcraft 

JMT-121  ^ . 1 

Subminiature  Switch,  Switchcraft 

JMT - 1 2 3 

Alpha  #3241  (2  conductor  w/shield) 
Alpha  #3242  (3  conductor  w/shield) 
Amphenol  Connector  7 Pm  #195 
(Cable  to  Panel  Plug) 

Power  Supply  Zeltex  ZM1550 
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' ....  . ....  ...^ 


:V’.  / '■ ‘■'1:..^.^....;'  ■'  ''' 


E„ut(T)  F.  G Rs(0) 


T-T 


ZRb+Rg 


a (1+6/100) 


(3) 


where  = Temperature  (°C)  when  bridge  is  nulled 


E = Bridge  bias  voltage  (-15  v) 
= Bias  resistance  (2.74  Kfi) 


R = Resistance  in  two  legs  of  the  bridge  (30.1  Kf2) 

a 


and  G ■ Voltage  gain  of  the  amplifier  (-5500) 


(In  evaluating  Eq.  (3),  we  approximated  Eq.  (2)  to  within  1%  for 
the  temperature  range  of  0 to  30°C.)  After  substitution  of  the 
component  values  into  Eq.  (3)  we  find 


E ^(T) 
out 


R,(0) 


T-T  108.2 

1 


(4) 


where  Rg(0)  is  in  ohms  and  ^^out ^ ^ v/®C.  Hence  a 

sensor  resistance  of  54.1  n will  yield  a 0.5  v/®C  gain  rf  the 
system.  (A  one  volt  per  ®C  gain  could  be  achieved  with  a value 
of  Rc(0)  « 54.1  n by  changing  R to  1 Mn.)  To  convert  an  actual 
output  voltage  to  temperature  changes,  multiply  by 


108.2/Rg(0) . 


Note  that  the  above  two  equations  have  Rg(0)  as  an  ele- 
ment. Normally  Rg  is  measured  at  some  temperature  other  than 
0®C.  So  it  is  convenient  that  soon  after  measurement  of  Rg,  it 
be  converted  to  Rg(0)  either  by  using  the  Calender-van  Dusen 
equation,  Eq.  (2)  or  the  numbers  in  Table  1. 
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The  maximum  operating  range  of  the  circuit  is  ± lOv, 
so  for  a gain  of  0.5  v/-C,  the  range  of  temperature  excursions 
is  . 20*C.  The  ambient  operating  point  can  drift  by  only  that 
amount  when  the  function  switch  is  in  the  DC  position.  However, 
when  AC  operation  is  used,  the  operating  point  may  drift  around 
in  the  range  i 200"C.  This  assumes  that  the  drift  is  of  fre 
quencies  below  the  1.6xl0->  Hz  high-pass  filter  cutoff.  Thus 
AC  operation  allows  the  sensor  to  be  run  untended  for  long 
periods  of  time.  It  also  allows  the  user  to  beep  the  tempera- 

A n V which  is  optimum  for  the  recording 
ture  fluctuations  around  0 v,  whicn  is  v 

and  subsequent  digitization  and  processing. 


C.  Operational  Testing. 


The  sensor  units  may  be  operationally  tested  in  two 
(We  will  assume  that  a well-calibrated  probe,  perhaps 
from  a commercial  concern,  is  not  available  as  a reference.) 

Both  approaches  involve  recording  output  data  and  processing  t e 

signals  for  power  spectral  densities.  The  processing  scheme 

first  test  is  for  the  entire 
used  is  outlined  in  section  V.  The  fir 

unit,  with  the  exception  of  the  probe  wire  itself.  The  pro  e is 
carefully  immersed  in  oil,  so  that  minimum  temperature  fluctua- 
tions are  sensed.  The  unit  is  operated  normally  otherwise. 

This  gives  us  a measure  of  the  noise  level.  Second,  we  operate 
a series  of  probes  in  the  atmosphere  under  normal  operating  con 
ditions.  Their  outputs  are  recorded  and  compared  with  theory 


and  each  other. 


The  oil  bath  test  consists  of  submersing  just  the  probe 
in  approximately  2 qt.  of  #10  weight  oil  and  allowing  the  tem- 
perature of  the  oil  to  stabilize  to  ambient  temperature.  With 
the  exception  of  the  probe  being  in  the  oil,  care  should  be 
taken  to  have  all  other  operating  conditions  as  nearly  normal  as 
possible.  That  is,  the  sun  should  be  shining  on  the  unit,  and 
the  wind  velocity  should  be  at  a reasonable  value.  For  an  addi- 
tional reference,  a second  probe  should  be  operated  normally  in 
the  atmosphere  within  a few  meters  of  the  sensor  being  tested. 
Once  confidence  has  been  developed  in  constructing  the  units, 
this  step  of  an  oil-bath  test  can  be  skipped.  An  example  of 
such  spectra  is  given  in  Fig.  4.  Total  mission  time  was  60  min 
and  Crp  - 0.113  m Note  that  this  value  of  turbulence  is 

high,  and  there  is  a two  decade  separation  between  the  two 
spectra  shown.  There  are  two  regions  of  the  oil  spectrum  to 
note.  The  lowest  frequency  power,  below  0.5  Hz  is  still  domina- 
ted by  drift  of  the  oil  bath;  hence  an  power  law.  The  re- 
mainder of  the  power  is  the  composite  noise  of  the  bridge,  am- 
plifier, coaxial  cable,  tape  recorder  amplifiers  and  analog 


tape.  The  conclusion  is  that  for  reasonably  high  values  of  C^, 


say  above  10  ^ m the  noise  is  not  a significant  prob- 

lem. If  we  include  the  drift  region  of  this  oil  spectrum,  the 
equivalent  noise  level  is  found  from  an  integral  of  the  spectrum 
to  be  0.016®C  (rms);  however  if  we  ignore  the  drift,  the  noise 
level  is  0.015®C  (rms).  The  latter  figure  is  consistent  with 
the  value  of  0.02®C  reported  by  Lawrence  et  al.® 
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The  second  test  to  be  performed  Is  to  look  at  the  out-  | 

puts  from  a series  of  probes  mounted  In  close  proximity  to  each  - 

other.  We  know  theoretically  what  the  temporal  spectrum  and 
spatial  covariance,  or  structure  function,  should  be.  Moreover 
the  variances  of  the  probes  should  be  within  IM  for  a one  hour 
mission.  It  is  convenient  to  mount  the  probes  on  a horltontal 
bar  and  to  space  them  in  a geometrical  progression.  That  pro- 
gression provides  the  largest  number  of  spatial  separations  be- 
tween the  various  probes.  This  tost  should  be  run  at  about  mid- 
day on  a day  with  mostly  clear  skies.  Height  above  ground  for 
the  horizontal  bar  should  be  a few  meters  so  that  local  ground 
effects  are  not  sensed.  All  separations  should  be  within  the 
outer  scale,  so  that  the  values  of  the  structure  function  fit 
the  Kolmogorov  form  given  by  E,.(l).  Exemplary  data  for  the 
structure  function  is  given  in  Fig.  5.  When  plotted  on  a log 
log  plot  as  shown,  the  inertial  subrange  should  appear  as  a 
straight  line,  which  when  extrapolated  to  the  unity  line  gives 
us  a value  of  the  outer  scale  . In  section  V of  this  report 
we  will  show  how  the  processing  for  this  structure  function  is 

done. 

Also,  for  each  sensor,  we  can  process  the  data  for  power 
spectra.  This  would  localize  problems  in  terms  of  frequency  and 
might  show  errors  ignored  in  the  structure  function  test.  For  | 

example  significant  60  Hz  ripple,  recorder  head  noise,  or  | 

amplifier  drift  might  be  isolated.  In  the  inertial  subrange,  | 

basically  between  v/(2vL^)  and  v/(2vt^).  the  spectrum  should  | 
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Normalized  temperature  structure  function  data 
collected  1 Oct  73  at  Verona  test  site,  Griffiss  AFB 
NY.  Outer  scale  L indicated  as  intersection  of  the 
extrapolation  of  tfie  inertial  subrange  with  the  unity 
line . 
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have  an  power  law  behavior  (except  perhaps  at  the  ends  of 

that  range).  An  example  of  such  spectra  corresponding  to  a 150 
minute  recording  time  and  two  probes  space  0.43  m is  shown  in 
Fig.  6.  Note  the  curve  is  quite  smooth  over  the  inertial  sub- 
range. How  we  obtain  such  spectra  is  also  outlined  in  the 
fiuai  section. 
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Figure  6.  Temperature  power  spectral  density  data  for  same 
mission  as  Fig.  5.  Two  probes  were  operated  simu 
taneously.  Inertial  subrange  indicated  as  an  £" ® 
power  law. 


' V- 


III.  SENSOR  PRORR  REPAIR,  CALIBRATION  AND  OPERATION 
A.  Manufacture  of  Sensor  Probes. 

The  fabrication  of  the  Wollaston  wire  probe  requires 
that  the  wire  be  soldered  to  a pair  of  sturdy  binding  posts  of 
brass  prior  to  etching  with  nitric  acid.  (Refer  again  to  the 
photograph  of  Pig.  2.)  The  binding  posts  are  themselves  sold- 
ered securely  to  a two  wire  cable  with  aluminum  shield.  The 
shield  is  soldered  only  to  the  printed  circuit  card  chassis 
ground  strap  to  minimize  60  llz  ground  loops.  Earlier  sensor 
units  were  fabricated  with  a probe  plug-socket  arrangement. 
However  an  investigation  of  Calender-van  Dusen  equation,  Eq.  (2), 
shows  that  only  a 0.004  rms  change  in  the  apparent  sensor  re- 
sistance duo  to  variations  in  the  probe-socket  resistance  will 
give  an  additional  0.02°C  rms  noise.  (This  assumes  R,.(0)  = 

54.1  n.)  Indeed  even  larger  variations  in  the  probe-socket 
resistance  were  noted  and  this  approach  was  abandoned.  This  em- 
phasizes the  need  for  good  solder  connections. 

At  this  point  we  assume  that  the  binding  posts  are 
available  and  will  describe  how  to  mount,  etch  and  test  the 
probe  wire.  Later,  we  will  describe  the  construction  of  the 

entire  unit.  The  following  steps  are  recommended  for  manufac- 
ture of  the  probes: 

1.  Prepare  the  nitric  acid  for  etching  in  this  manner: 

a.  Obtain  a 40?i  solution  of  reagent  grade  HNO  . 

3 

b.  Use  an  eye  dropper  to  produce  a controlled  drop  of  acid. 

20 
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4. 


wires  of  the  probe,  preferably  making  the  connections 
at  the  printed  circuit  card.  (An  HP  3440A  with  plug 
in  3444A  is  an  adequate  ohmmeter.) 

c.  Record  the  resistance  value  which  may  be  as  high  as 
0.4  n. 

d.  Place  about  2 mm  of  the  tip  of  the  wire  loop  in  the 
suspended  drop  of  nitric  acid  by  lowering  the  eye  drop- 
per. At  normal  room  temperature  it  will  take  about  ten 
minutes  for  the  probe  resistance  to  reach  approximately 
58  n.  At  lower  temperatures  (say  10®C),  etching  may  be 
extremely  slow,  and  an  auxiliary  heater  may  be  needed  to 
raise  the  ambient  temperature  to  about  20  to  25®C. 

e.  Remove  the  probe  wire  when  the  resistance  is  approxi- 
mately 58  (±  5 n typically).  Note  that  a 2 mm  length 

of  20  n/mm  wire  would  give  us  40  Q resistance;  so  we  are 
certain  that  the  silver  is  removed  if  we  go  beyond  40  n. 

f.  Remove  the  probe  support  rod  from  its  clamp  and  momen- 
tarily submerge  the  wire  loop  in  a breaker  containing  a 
saturable  solution  of  baking  soda  in  tap  water.  This’ 
neutralizes  any  residual  acid. 

Complete  the  calibration  of  the  sensor  with  these  steps: 

a.  Record  the  probe  resistance  and  the  ambient  temperature 
in  the  vicinity  of  the  probe. 

b.  Subtract  the  initial  resistance  value  as  found  in  step 
3.  c.  from  the  resistance  found  in  step  4.  a. 

c.  Use  Table  1 or  the  Calender-van  Dusen  equation  to 
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determine  given  R^(T),  where  T was  the  ambient 

temperature  at  time  of  measurement. 

Record  the  sensor  identification  number  and  the  value  of 
R^(0).  It  is  this  value  of  that  should  be  kept 

with  the  unit  until  the  sensor  breaks;  Rs(T)  is  not  of 
value  any  longer  and  may  confuse  those  using  the  data, 
connect  the  DC  supply  leads  to  the  unit,  set  the  func- 
tion switch  to  ground  then  to  the  DC  position,  and  ad- 
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just  the  bridge  balance  pot  R^  so  that  the  output  vol- 


£. 


tage  is  nearly  zero  (i  0.05  v) . A slight  imbalance 

here  will  not  affect  the  final  output,  as  we  will  oper- 
ate  in  an  AC  coupled  mode. 

Measure  and  record  the  rms  voltage  output  (with  the 
probe  covered).  This  value  should  not  exceed  2 mv. 


B.  Field  Operation. 


For  optimum  turbulence  response,  the  sensor  should  be 
oriented  so  that  there  are  no  obstructions  upstream  of  the 
probe.  If  the  wind  direction  is  fairly  steady  (10  rms),  this 
means  that  the  rod  on  which  the  probe  is  mounted  be  pointed  into 
the  wind.  If  the  wind  direction  fluctuates  considerably,  then 
the  probe  support  rod  should  be  mounted  vertically  to  have 
proper  exposure  for  all  wind  directions. 


If  an  array  of  probes  is  used,  so  as  to  develop  a spa- 
tial correlation  or  structure  function,  then  the  line  of  probes 
should  be  normal  to  the  wind  direction.  Then  there  will  be  no 
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wake  effects  from  one  sensor  to  the  next.  Generally  the  data 
will  still  be  acceptable  as  long  as  the  wind  direction  does  not 
come  within  25°  of  the  line  connecting  the  probes  for  more  than 
10^  of  the  time.  The  operator  should  be  prepared  for  large 
fluctuations  in  wind  direction,  and  if  the  mean  shifts  radi- 
cally, he  should  reposition  the  bar. 


The  distance  downwind  from  an  obstacle  such  as  a building 
should  be  at  least  100  times  the  size  of  the  obstacle.  This  in- 
cludes concrete  or  asphalt  roads  or  parking  lots  as  well.  Of 
course,  there  may  be  times  when  the  user  wishes  to  know  turbu- 
lence levels  in  the  vicinity  of  such  obstacles  and  this  rule  will 
not  apply.  Good  quality  microtemperature  data  will  be  collected, 
Iiowever,  when  well  - developed  turbulence  conditions  exist.  Typi- 
cally this  means  mostly  clear  skies,  daytime  conditions  and  wind 
speeds  exceeding  1 m/sec. 


The  sensor  power  supply  requirements  are  + 15  v and 
15  V at  ± 3 ma  maximum.  The  supply  may  be  a commercial  recti- 
fier with  filtering  or  may  b?^  a battery.  The  maximum  length  of 
the  power  cable  should  be  100  ft.  The  power  cable  is  of  #18 
wire,  including  2 twisted  leads  and  a ground  shield  lead.  Power 
supplies  may  be  located  near  the  sensors,  so  this  length  require- 
ment is  not  a real  restriction.  The  sensor  unit  signal  cable 
should  also  be  two  wires  plus  a ground  shield.  Normal  load  to 
the  sensor  is  10  Kfi.  a signal  cable  of  #18  wire,  one  mile  in 
length,  is  less  than  100  in  resistance,  lienee  the  signal  cable 
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Make  sure  there  is  no  voltage  saturation  occurring  in  the 
± 10  V operational  range. 

12.  The  sensor  unit  now  operational.  The  probe  should  not  be 
uncovered  in  winds  over  12  m/sec  or  in  conditions  oi'  dust  or 
precipitation,  as  they  will  break  easily.  In  fact,  any  sort  of 
jar  to  the  system  might  break  or  stretch  the  probe. 
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IV.  CONSTRUCTION 


In  previous  sections  we  considered  that  the  operational 
units  are  available  to  the  user  and  thus  we  only  covered  set-up, 
calibration  and  routine  usage.  Now  we  will  describe  how  the 
circuit,  enclosure  and  probe  support  are  manufactured,  especially 
emphasizing  rationale  for  certain  design  schemes.  In  retrospect 
we  feel  this  is  not  necessarily  a compact  design,  as  since  these 
units  were  first  built,  more  compact  integrated  circuitry  has 
been  made  available. 


A.  The  Circuit. 


Fabrication  of  the  high  gain  sensor  unit  requires  good 
engineering  practice  in  order  to  reduce  micx-ophonics , component 
value  drift,  rf  noise  pick-up,  60  Hz  pick-up  and  corrosion. 
Microphonics  have  been  eliminated  by  using  a printed  circuit 
board  and  strong  solder  connections.  Noise  pick-up  from  exter- 
nal rf  sources  is  greatly  reduced  by  mounting  the  printed  circuit 
card  inside  an  enclosed  aluminum  chassis,  and  by  mechanically 
attaching  the  probe  support  tube  to  the  chassis.  By  enclosing 
the  circuitry  in  the  aluminum  box,  no  air  currents  are  present 
which  could  induce  temperature  variation  in  tlie  components  and 
thus  component  value  drift.  Corrosion  of  components  and  solder 
joints  is  eliminated  by  using  weather  sealed  covers  on  the 
aluminum  chassis. 


The  60  Hz  pick-up  is  reduced  to  less  than  the  amplifier 
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background  noise  output  level  by  using  two  independent  ground 
connections.  The  two  systems  are  not  connected  at  any  point 
but  at  the  recording  system  so  as  to  minimize  ground  loops. 

(To  better  understand  the  following,  refer  to  P'ig.  7 of  the 
printed  circuit  card.)  The  first  ground  is  a connection  of  the 
shield  of  the  output  signal  cable  to  the  GROUND  RUSS  on  the 
printed  circuit  card.  This  in  turn  is  connected  to  the  aluminum 
chassis  and  the  shield  for  the  cable  leading  to  the  microthermal 
probe.  As  with  the  signal  cable,  a two- twisted-wire  cable  with 
shield  is  used  for  the  probe.  The  power  supply  cable  is  howevc-r 
three-wire  with  a shield,  but  this  shield  is  connected  at  the 
poA^'er  supply  end  of  the  cable  directly  to  the  power  supply  case. 
None  of  the  wires  from  the  power  supply  is  connected  to  the 
aforementioned  GROUND  BUSS.  Tlie  second  ground  system  is  a con- 
nection of  the  ''low''  sides  of  the  probe  and  output  signal  two- 
wire  loads  to  the  COMMON  GROUND.  The  power  supply  zero-voltage 
lead  is  also  connected  to  the  COMMON  GROUND.  Thus  to  avoid  a 
ground  loop  from  the  power  supply,  this  lead  must  be  floated 
at  the  power  supply. 


The  gain  of  the  sensor  unit  must  have  a high  order  of 
stability  over  a range  of  ambient  temperatures  from  -30®C  to 
40®C.  Therefore  all  resistors  are  metal  film,  of  ± 1%  toler- 
ance and  have  a temperature  coefficient  of  ± 100  PPM/“C  from 
-55“C  to  175°C. 


The  printed  circuit  board  layout  of  Fig.  7 is 
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Fig,  7a.  OBVERSE  VIE’.i 


SIGNAL 

CABLE; 


CAL  SW 

(Off  when 
shorted) 


PROBE 

CABLE: 

HI 

LOW 


SHIELD 


Fig.  7b.  REVERSE  VIEW,  where  the  Printed  Circuit  is  etched. 


Figure  7.  layout  of  printed  circuit  board.  Sym’^olr-  A,  B and  C 
denote  the  connection  points  for  the  function  switch. 
(Connecting  A to  B implies  GMD,  B to  C implies  DC,  and 
no  connection  for  AC.  See  Fig.  9.) 


w 
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supplemented  with  an  exact  1:1  scale  layout  of  the  mask  used  in 
photo-etching,  Fig.  8.  Note  that  the  printed  circuit  is  under- 
neath the  board  shown  in  Fig.  7,  and  the  components  are  on  top. 
Thus  the  mask  in  Fig.  8 is  reversed  from  the  layout  of  Fig.  7. 

In  Figures  9 and  10  are  the  obverse  and  reverse  side  views  of 
the  board  after  components  have  been  mounted,  and  after  the 
board  has  been  secured  to  the  chassis.  Note  that  components 
C , C , C and  R are  mounted  directly  to  the  operational  ampli- 
fier  sockets.  Five  small  additional  capacitors  show  up  in 
Figs.  9 and  10  and  were  added  just  prior  to  the  time  the  photo- 
graph was  taken.  They  are  not  listed  in  the  parts  list  of 
Table  2 and  are  not  shown  in  Fig.  7.  They  were  used  to  reduce 
exceptionally  high  rf  pick-up  in  one  experiment  and  connect 
from  the  GROUND  BUSS  to  ± 15  v,  0 v,  and  to  the  two  leads  of  the 
probe  cable.  Their  value  is  0.002  yfd  (50  WVDC) . 


B.  The  Mechanical  Support  Assembly. 


We  have  provided  full  engineering  drawings  of  the  sensor 
unit  in  Figs.  11-18  which  should  be  sufficient  for  the  construc- 
tion of  units  identical  to  ours.  Certain  points  may  be  improved 
upon;  however  our  intent  is  to  provide  plans  for  a working  system 
and  not  present  anything  not  tried  in  practice.  The  various 
photographic  views.  Fig.  1,  2,  9 and  10  will  be  helpful  in  fol- 
lowing this  discussion.  In  Table  3 we  have  listed  those  parts 
for  which  government  stock  numbers  exist.  Note  that  this  list 
does  not  cover  all  parts  needed,  as  we  anticipate  that  the 
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tographic  viev;  of  printed  circuit, 
hassis . 


components 


ures  11-18 


Engineering  drawings  of  sensor  unit,  includi 


Covers  for  chassis  ' 
Cap  to  cover  probe 
Sleeve  to  cover  probe 
Rods  for  mounting  probe 
Holder  for  rods  of  Fin. 


support  tube 
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Table  3. 


List  of  materials  for  constructing  sensor  unit, 
excluding  electronics. 


1.  Muminutn  tubing  - Stock  M7UI-278-6411,  1/2  inch  outside 
diameter,  wall  thickness  0.05  inch,  weldable  6061. 

2.  Aluminum  bar  - Stock  //9530-247-5514,  1/8  inch  thick,  2 inch 
wide,  weldahle  6061. 

3.  Aluminum  har  - Stock  //9530-237-0721  , 1/4  inch  thick,  1 1/4 
inch  wide,  weldable  60'ol. 


4. 


5. 


Aluminum  bar  - Stock  /I9530-244-9026  , 3/4  inch  thick,  3/4 
inch  wide,  weldahle  6061.  (Alternative:  aluminum  tubing 

3/4  inch  OH,  1/8  inch  wall  tliickness.) 


Aluminum  sheet  - Stock  /t‘9535-232-6868 , 0.063  inch  thick, 
weldable  5052. 


6.  Rubatex  Insulation  Tape  - 1/8  inch  thick  by  2 inch  wide 
(Rubatex  Corp.,  Bedford  VA) . 

7.  Teflon  (or  Nylon)  cylindrical  bar  - 1 inch  diameter 

8.  Teflon  (or  Nylon)  cylindrical  l)ar  - 1/2  inch  diameter 


i 

I 

I 

machinist  will  have  common  materials,  such  as  screws,  nuts,  etc.,  I 
in  stock.  5 

Rather  tlian  discuss  oacli  of  tlie  ciglit  engineering  draw- 
ings in  detail,  we  will  break  tlie  unit  into  three  general  ])arts: 
chassis  (Figs.  12,  13),  probe  cover  (Figs.  14,  15)  and  probe  i 

support  tube  (Figs.  16-18).  Tlie  avcrall  unit  is  assembled  in  I 

Fig.  11.  ^ 

The  cliassis  is  fabricated  fiom  weldable  aluminum  6061  . ; 

Chassis  sides  (long  dimension)  are  from  1/8  inch  by  2 inch  bar  i 

stock  and  the  ends  arc  from  1/4  inch  by  1 1/4  inch  bar  stock.  1 

The  chassis  corners  are  secured  by  welding,  The  two  chassis  i 

covers  are  fabricated  from  aluminum  flat  sheet  0.06  inch  thick. 

Edges  have  been  bent  down  for  a better  appearance  and  to  avoid  1 

sharp  edges.  Along  the  inner  edges  of  the  covers,  where  the  1 

covers  mate  to  the  chassis  sides,  strips  of  rubber  have  been  | 

used  so  that  air  currents  may  not  enter  the  box.  The  probe 
tube  socket  (labeled  in  Fig.  9)  may  be  made  from  3/4  inch  by 
3/4  inch  bar  stock,  and  is  welded  to  the  chassis.  An  alterna-  ! 

tive  material,  shown  in  Fig.  12,  is  stock  tubing,  0.75  inch  | 

outside  diameter,  0.5^  inch  inside  diameter.  Tlie  probe  tube  | 

socket  is  drilled  and  tapped  on  two  sides  for  #20  set  screws 

which  secure  the  tube.  The  tube  could  be  welded  in  place  if  the 

1 

user  never  expected  to  interchange  Iciger  or  shorter  tubes,  and  | 

1 

this  would  prevent  moisture  from  entering  the  chassis  at  this 
point.  The  tube  and  socket  could  be  threaded,  as  well.  The 


I 


printed  circuit  card  is  fastened  to  the  inside  of  the  cliassis 
by  means  of  two  cubical  aluminum  liolders  which  arc  tapped  for 


6/32  incii  screws.  The  holders  are  secured  to  tlic  inner  \v'alls 
of  the  chassis  by  welding,  as  shown,  h'c  suggest  preparing  tem- 
plates for  use  in  drilling  the  lioles  into  t!ie  chassis  for  the 
mounting  of  sv.'itclics  and  sockets,  and  for  tlic  holes  for  tliC 
chassis  cover  plates. 

Tlic  probe  support  tulic  (indicated  in  fig.  f and  drawn  in 
Fig.  18)  is  also  from  weldable  aluminum  stock  6061,  0 . f)  inch 
outside  diameter  with  a wall  thickness  of  0.049  inch.  At  the 
probe  end  of  the  tube,  a slot  is  milled  longitudinally  (along 
the  rod)  beginning  at  a point  1/2  inch  from  the  end  of  tne  tube 
and  extending  for  2 1/4  inches.  The  groove  is  J/8  inch  wide  and 
0.03  incli  deep.  In  addition  a 1/2  inch  long  slot,  1/8  inch  in 
width  and  0.03  inch  deep  is  milled  laterally  from  the  end  of  the 
longitudinal  slot  at  tlic  end  of  the  tube  nearest  the  probe. 

This  groove  ])rovidcs  a track  for  tlie  set  screw  in  the  teflon 
cover,  when  tlie  cover  is  liroughc  forward  to  cover  the  probe. 

The  lateral  slot  allows  the  cover  to  be  rotated  about  90°  after 
it  is  slid  over  the  ]irolic.  Tlie  sensor  end  of  the  probe  support 
tube  is  fitted  with  a teflon  cylinder  to  hold  the  two  brass  rods 
which  in  turn  support  the  Wollaston  process  wire.  (This  cylinder 
is  not  shown  in  the  photographs  but  is  detailed  in  Fig.  17.) 

The  cylinder  is  of  diameter  0.4  inch  and  is  1/2  inch  in  length. 

To  secure  the  cylinder,  a #4  set  screw  is  inserted  through  the 
side  of  the  tube.  The  rods  are  two  inches  long  and  75  mil  in 
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diameter  (re  PiR.  16).  The  teflon  cylinder  is  drilled  for  the 
Insertion  of  the  rods  at  the  spacing  shown  in  '-ig.  17.  Once  the 
rods  are  in  place,  the  prohe  cables  may  be  soldered  to  the  rods 
and  the  cylinder  inserted  into  the  tube.  The  material  used  for 
the  rods  is  a brass-alloy  welding  rod. 


The  probe  cover  is  fabricated  from  1 inch  outside  di- 
ameter solid  teflon  or  nylon  rod.  The  three  inch  long  section 
of  teflon  rod  shown  in  Fig.  15  is  drilled  for  smooth  insertion 
of  the  1/2  inch  diameter  tube.  The  rod  is  cut  back  0.38  inch 
of  its  length  to  an  outside  d:.ameter  of  3/4  inch.  A 1 inch 
diameter  by  1/2  inch  long  cap  (in  Fig.  14)  is  then  machined  to 
fit  over  the  3/4  inch  diameter  end  of  the  teflon  rod.  The 
opposite  end  of  the  rod  is  tapped  for  a #6  thumb  screw.  This 
screw  permits  the  probe  to  be  covered  by  sliding  the  teflon 
cover  over  the  end  of  the  sensor  and  attaching  the  cap.  The 
cap  is  fitted  with  a long  string,  as  shown  in  Fig.  1,  so  that 

the  cap  may  be  held  in  place. 
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V.  DATA  PROCESSING  SYSTEM 


A.  The  Recording  Scheme. 

The  purpose  o£  the  data  processing  system  is  to  provide 
temporal  and  spatial  statistical  information  on  the  turbulent 
temperature  of  the  atmosphere.  In  addition  to  the  microtempera- 
ture measurements,  common  meteorological  observables  which  are 
recorded  include  wind  speed,  wind  direction,  barometric  pressure 
and  the  gross  temperature  from  a dew-point  system.  Nine 
microthermal  probes  are  oriented  in  a geometrical  progression 
on  a horizontal  boom  positionable  two  to  eight  meters  above  the 
ground.  This  makes  a total  of  thirteen  signal;  which  are 
recorded  digitally  on  a narrow  bandwidth  system  and  on  analog 
tape  for  subsequent  wide  bandwidth  processing.  The  digital 
unit  is  the  Analog-Digital  Data  System  model  lOOM  Data  Logger, 
with  a 12  bit  analog-to-digital  converter.  The  sample  rate 
is  400  Hz  divided  by  the  number  of  channels  currently  in  use. 

For  the  13  channels,  the  rate  is  about  31  Hz.  This  narrow 
bandwidth  system  is  used  for  obtaining  those  statistical 
quantities  which  are  not  sensitive  to  aliasing:  especially 
probability  densities,  variances,  spatial  cross -covariances 
and  higher  order  moments.  For  spectral  information  on  the 
turbulence,  all  thirteen  signals  are  recorded  in  parallel 
on  analog  tape.  The  Ampex  FR1200  tape  recorder  has  at 
3 3/4  ips  a bandwidth  of  DC  to  1250  Hz  (-3  dB) 

and  the  tapes  used  have  bandwidths  superior  to  that.  Later  two 
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channels  at  a time  are  converted  from  analog-to-digital  on  a 
12-bit  PDP-8  computer.  Because  the  possibility  of  aliasing  will 
be  significant  to  the  subsequent  Fourier  analysis,  a low-pass 
filter  is  inserted  in  each  channel  prior  to  conversion.  The 
-3  dB  point  of  the  Krohn-llite  model  3322  filters  is  set  to  the 
Nyquist  frequency,  one-half  tlie  sample  rate,  and  the  roll-off 
is  that  of  a four-pole  Butterworth  filter,  -80  dB/decade. 

Sample  rates  from  200  to  700  llz  are  chosen  depending  on  the  ex- 
pected noise  level. 


B.  Digital  Processing 


The  digital  tapes  from  the  ADDS  unit  and  from  the  PDP-8 
are  processed  in  tv^ro  separate  software  packages  on  a Honeywell 
635/645  compu.ter.  One  program  computes  the  cross-covariance 
matrix  of  the  thirteen  signals  by  operating  on  the  tapes  from 
the  ADDS  unit.  Matrices  are  generated  for  averaging  times 
typically  ranging  from  10  to  30  minutes.  The  microtemperature 
elements  of  the  matrix  are  then  converted  into  a plot  of  1-C(.r) 
where  C(r)  is  the  spatial  correlation.  Such  a plot  is  identical 
to  D.j,(r) /D.j^  («>)  , v\rhere,  once  again,  Dx  is  the  structure  function. 
In  the  end  an  average  l-C(r)  plot  is  generated  for  the  entire 
mission.  This  program  also  calculates  C|  according  to  its 
definition  Eq.  (1)  and  calculates  as  the  separation,  r where 
the  extrapolated  inertial  subrange  of  l-C(r)  intersects  the 
unity  line.  An  example  of  such  data  is  given  in  Fig.  5. 


The  second  program  operates  on  the  tapes  from  the  PDP-8 
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computer  and  provides  wide  bandividth  power  spectral  densities. 


•or  a full  description  of  this  program  we  direct  the  reader  to 


Jayne  and  Foley®.  The  six  spectra  are  the  individual  spectra 


of  the  two  temperatures,  the  difference  spectrum,  the  modulus 
and  phase  of  tlie  cross -spectrum  and  the  coherency.  The  signals 


are  first  transformed  via  a fast  discrete  Fourier  algorithm  for 


frequencies  from  f^y/2048  to  f^^ , where  f^^  is  the  Nyquist  fre- 


quency. Namely  4096  data  points  per  chtinnel  are  transformed  at 
one  time.  The  data  are  then  digitally  low-pass  filtered  by  a 


four-pole  Butterworth  with  a cut-off  of  f^^/C4.  The  low-fre- 


quency part  of  the  spectrum  is  tlien  obtained  by  transforming  the 
remaining  signal  which  lias  been  decimated,  keeping  only  one  point 
in  sixty- four.  The  low  and  high  frequency  spectra  undergo  con- 
stant percentage  bandwidth  averaging  so  that  on  a logarithmic 
frequency  scale,  tlie  spectral  estimates  appear  to  have  a nearly 
constant  separation.  Sixty-four  frequency  spectra  are  averaged 
to  correspond,  in  averaging  time,  to  a single  low-frequency 
spectrum.  Finally,  several  such  spectra  are  averaged  to  repre- 
sent an  entire  mission.  The  final  averaging  time  can  b'  computed 


as  N*2048*64/fj^y,  wliere  N is  the  number  of  spectra  averaged  to 


gether.  Since  varying  numbers  of  spectral  estimates  are  aver- 
aged to  give  a final  value,  the  confidence  limits  are  a function 


of  frequency  as  well  as  N.  Following  Blackman  and  Tukey’,  we 
assume  the  spectral  estimates  follow  a chi-squared  distribution. 
Tables  in  Blackman  and  Tukey  or  in  Abramowitz  and  Stegun®  can 
be  used  to  generate  upper  and  lower  confidence  limits.  The 


confidence  limits  for  our  spectra  never  exceeded  twice  the  value 
at  the  high  side  or  1/2  the  value  at  the  low  side  and  the  error 
bar  size  decreased  rapidly  with  increasing  frequency. 

C.  Pertinent  Equations. 

We  now  assume  that  power  spectral  densities  and  covari- 
ances have  been  computed  and  will  demonstrate  the  important 
parameters  which  describe  turbulence.  Most  important  of  these 
parameters  are  the  temperature  structure  parameter  c“  and  the 
outer  scale  L..  We  will  show  how  to  insure  the  nature  of  the 
turbulence  is  Kolmogorov  and  how  to  treat  and  Lo  i£  it  is 

not . 

There  would  be  no  need  to  measure  both  the  temporal 
spectra  and  the  spatial  covariance  if  the  assumption  of  frozen 
flow  were  valid.  Frozen-flow  assumes  that  we  can  convert 
spatial  lags  into  temporal  simply  by  considering  the  wind 
velocity  as  a scale  factor.  If  the  turbulence  does  not  change 
significantly  while  it  is  transported  through  the  spatial  lags 
of  Interest,  then  the  frozen  flow  assumption  is  valid.  Panofsky. 
Cramer  and  Rao’  state  the  hypothesis  is  valid  close  to  the 
ground  for  lag  distances  up  to  90  m provided  i 1/J.  where 

0/  is  the  wind  speed  variance.  This  translates  to  equivalency 
of  temporal  and  spatial  spectra  for  temporal  frequencies 

£ ? v/570. 

So  why  would  we  measure  both  the  spatial  covariance  and 
the  temporal  power  spectra?  First  of  all,  we  are  not  guaranteed 
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frozen  flow  behavior.  Secondly,  we  can  not  obtain  an  accurate 

measure  of  outer  scale  from  the  spectra.  Thirdly,  to  obtain  a 
2 

value  of  C.p  from  the  spectra  requires  knowledge  of  wind  velocity; 
whereas  it  can  be  derived  early  from  the  covariance  of  two  probes 
with  only  the  knowledge  of  their  separation.  We  usually  use  only 
two  probes  and  calculate  their  spectra  and  a single  point  on  the 
structure  function. 

Let  us  first  examine  the  structure  function.  We  compute 
the  covariance  matrix 


CS) 


for  an  array  of  probes  numbered  i (or  j)  = 1,  N.  The  on-axis 
components  are  the  variances 


a • = C • • . 
1 11 


(6) 


For  the  structure  function  we  return  to  Eq.  Cl) » repeated  here 
for  convenience: 


n 


= <ni-Tpn 


(7) 


where  r^.  is  the  separation  of  the  two  sensors  at  i and  j.  If 
we  expand  this  equation  and  substitute  Eqs.  C5)  and  C6)  we  find 


ILpCr)  = a?  + aj  - 


C«) 


where  y^  = <T^>. 

Without  any  loss  in  generality  we  can  assume  • In  practice 

the  mean  values  may  be  offset  only  by  the  electronics,  as  long 
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fcr  ; 
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as  the  two  probes  are  at  the  same  height  above  ground.  Thus 


D (r.  .) 
T 13 


+ 0 - 2C . . . 

i 3 13 


C9) 


Z. 

We  are  now  in  a position  to  get  values  o£  C.^,,  as  long  as 


(by  typically  a factor  of  10).  In  a moment  we  will  see  how  to 


obtain  Lq . So  from  Eq.  (1)  we  have  as  many  measures  of  as  we 


have  r 


13 
-.2  _ 


f2  _ n fr--3r'?'^^.  r--<<l  . 

''T  iJ  ’ 0 


(10) 


We  do  not  have  to  concern  ourselves  with  r^^  being  so  small  as 
to  provide  a measure  of  the  inner  scale,  since  sensors  should 


not  be  placed  within  a few  millimeters  o£  each  other:  each 

would  interfere  with  the  flow  seen  by  the  other. 

If  the  sampling  is  of  sufficient  duration  (at  least  60 
minutes)  and  the  probe  units  are  accurately  calibrated  with 
respect  to  each  other,  then 


q = aj. 


(11) 


(In  addition  the  probes  must  be  at  the  same  height  above  ground.) 


Next,  we  assume  that  C..(“)-0,  which  is  true  for  all  physical 

13 


processes.  Then  if  we  divide  D.j.(r..)  by  D.j.(oo)  we  obtain 


b.p(rij)/IXj.(”)  = 1 - Pij 


(12) 


where  p.  is  the  cross-correlation  coefficient  of  the  two  sen- 

13 


sors.  It  is  the  function  I'Pj^j  which  is  plotted  in  Fig.  5 and 
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which  provides  us  with  the  outer  scale.  By  definition’*'  the 
outer  scale  is  the  separation  r at  the  intersection  of  an 


extrapolation  of  the  r^/^  region  and  the  1-p  = unity  line 

ij 


indicated  in  Fig,  5j  and  there  we  found  ~ 6 m for 
probes  2.2  m above  ground  level.  In  order  to  obtain  the  number 
of  separations  shown  in  Fig.  5,  we  oriented  nine  probes  in  a 


geometrical  progression.  If  we  have  m probes  in  such  a progres- 

n-1 

sion,  then  we  have  Z (n-i)  independent  separations. 

i = l 


The  other  piece  of  information  we  can  obtain  is  whether 
the  turbulence  is  Kolmogorov  or  not.  The  power  law  dependency 
in  the  inertial  subrange  must  be  r / . Our  rule  of  thumb  is  that 
the  least  squares  power  law  should  be  well  within  101  of  2/3 
for  recording  times  of  at  least  one  hour.  Also,  we  only  con~ 
sider  separations  r<!Lj/10  in  making  this  judgement.  The  values 


of  1-p..  for  small  r. . have  less  confidence  than  those  of  larger 
^ij  ' reason  is  that  if  we  assume  a fixed  percentage  confi- 

dence limit  on  p.  . , then  as  p.  .-»-l,  the  values  of  1-p  become 

ij  11  ij 

less  trustworthy.  If  we  determine  that  the  nature  of  the  turbu- 


lence is  non-Kolmogorov,  we  do  not  throw  out  our  value  of  L^. 

It  still  represents  an  "external- scale  size."  However  would 
be  recomputed  for  a different  structure  function  power  law  r^ 


and  would  be  given  a new  symbol  O'* 


The  relation  for  C*  is 

T 


p 

The  outer  scale  used  here  is  equal  to  the  one  obtained  from 
the  von  Karman  spectrum  if  you  multiply  the  von  Karman  L„  bv 
1.071.  ° 
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r*2  = n (r)r'l\  r<<L*, 

T '1'  ^ 0 


(13) 


v^here  L*  is  the  new  external  scale  size.  Note  that  i£  only  two 


probes  are  used  and  the  outer  scale  is  sufficiently  large,  then 
it  is  convenient  to  set  r = 1 m in  advance. 

The  power  spectra  also  demonstrate  whether  the  turbulence 
is  Kolmogorov.  Assuming  frozen  flow  and  the  form  of  the  struc- 
ture function  given  in  Eq.  (1),  we  find  the  spectrum  in  the  in- 
ertial subrange  to  be 


F^(f)  = 0.0730V  / f-V^C 


(14) 


for 

f * 

V 

2 2 
<<f  << 

V 

2ttLo 

2ttJIo 

1 

, where  v = wind  speed.  Note  that  we 


now  have  another  way  of  measuring  and  it  is  from  a single 


probe.  However  we  now  have  the  added  requirement  of  measuring 
wind  speed.  For  convenience  we  rewrite  Eq.  (14): 


= 13.69  f5/3v-2/3F^(f) 


(15) 


Values  of  outer  scale  could  be  obtained  from  spectra  but  it  re- 
quires recording  times  typically  much  longer  than  periods  of 
stationarity  of  the  turbulence.  It  is  much  more  convenient  to 
use  a pair  of  probes  at  a fixed  separation  and  measure  their 


correlation  coefficient.  An  example  of  obtaining  from  spectra 


is  shown  in  Fig.  6.  The  inertial  subrange  power  law  dependence 

is  shown,  and  from  a value  of  F^  at  •£=!  -Hz  and  v=3.2  m/sec 

2 
T 


we  find  Ct  = 0.20  °C^m-V\ 
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ihe  power  spectra  can  be  used  to  give  measures  of  inner 
scale  Zg.  However  the  sensors  used  must  be  significantly  shorter 
than  Zg,  so  as  not  to  effect  a spatial  average.  Furthermore, 
the  sensor  must  respond  well  in  frequencies  around  v/(2nZg)  which 
may  be  a few  hundred  Hz.  This  may  require  a very  fine  wire  sen- 
sor of  diameter  0.6  ym.  Our  sensors  respond  linearly  to  approxi- 
mately 100  Hz,  are  of  diameter  2.5  ym  and  are  2 mm  in  length. 
Hence  they  are  not  quite  suitable  for  measuring  Zg.  The  only 
reason  we  do  not  use  finer,  shorter  probe  wires  is  that  they  are 
difficult  to  construct  and  break  easily. 

The  quantity  of  interest  in  the  propagation  calculations 
IS  C^,  the  refractive- index  structure  parameter.  Once  we  make 
the  assumption  that  water  vapor  and  pressure  fluctuations  are 
insignificant,  we  find  that  is  proportional  to  via 


= C79.2x10-®p/t')'  c' 


C16) 


where  p is  the  ambient  pressure  in  millibars  and  T is  the  am- 
bient temperature  in  degrees  Kelvin.  The  constant  79.2x10'®  is 

correct  only  at  X=0.5  ym,  in  the  middle  of  the  visible  range. 

1 0 

Tatarsldi  also  reports  these  values  for  the  constant; 

82.9x10-'  at  1=0.3  ym  and  77.4x10-'  at  X=0.7  ym.  We  have  per- 
formed a least  squares  fit  to  these  three  values  of  the  constant 
and  found  the  relation 


2 

(72.3503  + 0.6067k  - 0.004912k  ) x 10 


- 5 


(17) 

where  k=2;r/X  and  X = wavelength  {um"‘)  in  the  visible  range.  We 
suppose  the  relation  is  not  too  far  off  in  the  infrared  but  have 
no  values  to  compare  with  in  that  range. 
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VI.  CONCLUSIONS 


We  have  presented  the  details  necessary  to  constructing 
and  using  a probe  for  measuring  microtemperature  fluctuations 
in  the  atmosphere.  Tlie  design  is  for  an  inexpensive  unit  which 
may  cost  only  about  $200  as  compared  with  commerical  ones  at 
least  ten  times  that.  Our  circuit  has  no  automatic  compensation 
should  the  level  of  fluctuations  change  radically,  and  there  is 
no  automatic  gain  circuit  which  would  compensate  for  low  fre- 
quency drifts.  The  unit  is  expected  to  be  tended  to  the  extent 
that  it  is  to  work  at  near  ground  level  and  may  be  covered  in 
case  of  precipitation. 


Also  we  have  included  no  internal  processing  of  the 
signal.  The  data  is  raw  temperature  fluctuations,  wliich  is  pre- 
ferred in  the  type  of  processing  we  do.  Should  the  user  decide 
to  redesign  the  unit,  he  should  keep  in  mind  the  major  design 
criteria  outlined  throughout  the  text.  V/e  expect  that  with  more 
up-to-date  circuitry,  the  entire  bridge-amplifier  arrangement 
might  be  placed  in  a volume  the  size  of  the  teflon  probe  cover. 
If  we  were  marketing  this  device,  that  is  what  we  might  have 
done.  However,  we  have  reported  on  a reliable  instrument  which 
has  served  us  for  six  years,  and  we  did  not  change  anything  in 
tlie  design  just  for  the  purpose  of  this  report. 
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